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A multi-phase, experimental study in the Basic Aerodynamics Research Tunnel at the 
NASA Langley Research Center has provided new insight into the unsteady flow interaction 
around cylinders in tandem arrangement. Phase 1 of the study characterized the mean and 
unsteady near-field flow around two cylinders of equal diameter using 2-D Particle Image 
Velocimetry (PIV) and hot-wire anemometry. These measurements were performed at a 
Reynolds number of 1.66 x 10 5 , based on cylinder diameter, and spacing-to-diameter ratios, 
L/D, of 1.435 and 3.7. The current phase, Phase 2, augments this dataset by characterizing 
the surface flow on the same configurations using steady and unsteady pressure 
measurements and surface flow visualization. Transition strips were applied to the front 
cylinder during both phases to produce a turbulent boundary layer upstream of the flow 
separation. For these flow conditions and L/D ratios, surface pressures on both the front 
and rear cylinders show the effects of L/D on flow symmetry, pressure recovery, and the 
location of flow separation and attachment. Mean streamlines and instantaneous vorticity 
obtained from the PIV data are used to explain the flow structure in the gap and near-wake 
regions and its relationship to the unsteady surface pressures. The combination of off-body 
and surface measurements provides a comprehensive dataset to develop and validate 
computational techniques for predicting the unsteady flow field at higher Reynolds 
numbers. 


Nomenclature 

= cylinder diameter 
= drag coefficient, per unit span 

= pressure coefficient, — — — 
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(7 p 

= fluctuating pressure coefficient, 
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= distance between cylinder centers 
= pressure measured on model surface 
= freestream static pressure 
= freestream dynamic pressure 

= Strouhal Number, 

U x 

= freestream velocity 
= streamwise coordinate 
= vertical coordinate 
= spanwise coordinate 
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AZ = distance between transducer and reference transducer for spanwise pressure correlations 

6 = azimuthal position 

o P = pressure standard deviation 


dV dU 

co z = spanwise vorticity component 


I. Introduction 

For the past several years, researchers at NASA have been performing experimental and computational studies 
on tandem cylinder arrangements to better understand the noise generation process and develop accurate noise 
prediction methods for landing gear systems. The tandem cylinder arrangement produces the type of unsteady flow 
interaction commonly associated with landing gear systems where upstream components, directly exposed to the 
oncoming flow, act as bluff bodies and generate unsteady wakes that impinge on one or more downstream 
components. As such, it provides a canonical problem to advance modeling techniques for the prediction and 
analysis of unsteady flow interactions. 

Numerous studies have been conducted on the tandem configuration due to its practical application in systems 
used for cooling, venting, and structural support. In most cases, the studies were performed at subcritical Reynolds 
numbers (less than 1.5 x 10 5 ) and focused on fundamental issues such as boundary layer development, numerical 
prediction of flow interference, flow-induced vibrations, and characterization of the wake structure for various 
cylinder spacings. 1 " 4 In Reference 5, Zdravkovich characterizes flow patterns associated with various separation 
distances (L/D) between the cylinders. For L/D < 1.1, the two cylinders behave as a single bluff body with vortex 
shedding occurring at the rear cylinder only. For 1.1 < L/D <1.6, the shear layers from the front cylinder attach 
alternately to the rear cylinder and vortex shedding occurs on the rear cylinder. For 1.6 < L/D < 2.5, the shear layers 
from the front cylinder attach to the rear cylinder and vortex shedding occurs on the rear cylinder. For separation 
distances, 2.5 < L/D < 3.2, intermittent shedding can be detected in the region between the cylinders and vortex 
shedding occurs on the rear cylinder. For 3.2 < L/D < 3.8, the flow between the cylinders is bistable and switches 
between intermittent shedding and constant shedding. For L/D >3.8, vortex shedding occurs on both cylinders with 
the same characteristics as a single cylinder. Quantitative details about the time-averaged and instantaneous flow 
structure for several of the flow regimes described above have been obtained by Lin, et al. 6 at Re = 1.0 x 10 4 using 
Particle Image Velocimetry (PIV). 

For supercritical Reynolds numbers, Gu 7 has shown from static pressure measurements that two flow regimes 
can exist for tandem cylinders depending on whether the rear cylinder is located in or out of the wake cavity of the 
front cylinder. For L/D < 2.8, the rear cylinder is located in the wake cavity and the shear layers from the front 
cylinder either envelope or impinge directly on the rear cylinder. For L/D > 2.8, the rear cylinder is out of the wake 
cavity and the shear layers from the front cylinder do not impinge directly on the rear cylinder. These descriptions 
provide useful information about the time-average flow between the cylinders but quantitative information about the 
unsteady flow field is also needed to enable computer model assessment and validation at critical and post-critical 
Reynolds numbers. Most studies of the unsteady flow field between tandem cylinders have used fluctuating surface 
pressure measurements and flow field spectra obtained using hot-wire anemometry to examine the effects of 
Reynolds number, freestream turbulence, surface roughness, and cylinder spacing on the flow structure associated 
with the different flow regimes described above. 8 " 11 Recent studies by Jenkins, et al. 12 , Khorrami, et al. 13 , and 
Henderson and Brooks 14 have added to this body of knowledge by providing additional insight into the unsteady 
near-field flow structure and its implications on noise radiation. 

The work by Jenkins, et al. 12 and Khorrami, et al. 13 constitute Phase 1 of a multi-phase study to measure and 
characterize the unsteady flow on two cylinders in tandem arrangement at a post-critical Reynolds number. The 
purpose of the current phase, Phase 2, is to document the steady and unsteady surface flow. The results from the 
current phase, along with those from Phase 1, provide a comprehensive dataset for 1) characterizing the interaction 
of the wake of the front cylinder with the rear cylinder at two different L/D’s, 2) examining differences between the 
wake of the entire tandem configuration and a single cylinder due to this interaction and 3) developing and 
validating computational techniques for the prediction of unsteady flows. 
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II. Experimental Apparatus and Techniques 


A. Test Facility 

The experiment was conducted in the NASA-Langley Basic Aerodynamics Research Tunnel (BART). This 
tunnel is a subsonic, atmospheric wind tunnel used to investigate the fundamental characteristics of complex flow 
fields and to acquire detailed data for the development and validation of computational fluid dynamics (CFD) 
models and methods. The tunnel has a closed test section with a height of 28 inches (0.71 1 m), a width of 40 inches 
(1.016 m), and a length of 120 inches (3.048 m). The free stream velocity was set to 184 ft/s (56 m/s) during Phase 
1 and 144 ft/s (44 m/s) during Phase 2 to achieve a Reynolds number based on cylinder diameter of 1.66 x 10 5 . At 
these conditions the free stream turbulence level is less than 0.10%. Additional information about the BART can be 
found in References 15-17. 

B. Model 

The tandem arrangement used for both phases of this study consisted of two circular cylinders which spanned the 
entire height of the test section. Figure 1 shows an installation photo and Fig. 2 provides a schematic of the models 
with the coordinate reference system. For Phase 1, each cylinder measured 1.75 inches (44.45 mm) in diameter and 
27.188 inches (690.6 mm) in length and was attached to a base that was inserted into a track mechanism secured to 
the tunnel floor. With this arrangement, the front cylinder was fixed and the rear cylinder could be translated 
upstream and downstream to achieve the desired spacing. In addition, data on a single cylinder could be acquired by 
fixing the location of the front cylinder and removing the rear cylinder. Each cylinder was instrumented with lower, 
center, and upper static pressure rings located at 10 inches, 14 inches, and 18 inches (254.0 mm, 355.6 mm, and 
457.2 mm), respectively, from the tunnel floor. The upper and lower rings consisted of 18 orifices spaced every 20 
degrees around the cylinder circumference, whereas the center ring consists of 36 orifices spaced every 10 degrees 
around the cylinder circumference. In addition, 9 pressure orifices with 2 inch (50.8 mm) spacing were installed 
along the span of each cylinder at 0= 180° to monitor pressure recovery and check for spanwise variations. 

For Phase 2, the same arrangement of static pressure orifices was used but in addition, 22 piezoresistive, 
differential pressure transducers were installed in each cylinder. To accommodate the pressure tubing and electrical 
wiring, two new cylinders measuring 2.25 inches (57.15 mm) in diameter were manufactured for use with the same 
track mechanism used in Phase 1. Increasing the diameter of the model increased the tunnel blockage to 5.63% 
compared to 4.37%, the blockage for Phase 1. The transducers were mounted inside the cylinders in two arrays: a 
spanwise array of 15 transducers along the length of the cylinder at 0-135° and a circumferential array of eight 
transducers at the center of the cylinder evenly spaced at 45° and including the transducer from the longitudinal 
array at 0 = 135°. To minimize flow disturbance, each transducer was mounted beneath the surface of the cylinder 
under a pinhole of approximately 0.05 inches (1.270 mm) in length and 0.02 inches (0.508 mm) in diameter. The 
fabrication process for installing the transducers in each cylinder left a conical chamber between the transducer face 
and the pinhole leading to the surface. This geometrical configuration raised concern about data contamination due 
to a Helmholtz-type resonance 18 . Calculations indicated that this type of resonance would likely occur around 14 
kHz and subsequent checks on a sample configuration in the NASA Langley Flow Impedance Test Facility 
confirmed that the resonance frequency was higher than the frequency range of interest for this experiment. 

To produce a turbulent boundary layer and achieve the type of flow separation associated with post-critical 
Reynolds numbers, transition strips were attached along the entire span of the upstream cylinder in the regions 
between 0 = 50° and 0 = 60° and between 0 = 300° and 0 = 310°. The transition strip was made of “pinked” tape 
with a nominal thickness of 0.229 mm. “Pinked” denotes that one edge of the tape was “straight” while the other had 
a “saw-tooth” pattern. The strip was applied with the “straight” edge facing upstream and the “pinked” edge facing 
downstream. In Fig. 3, the effectiveness of the transition strip in producing a post-critical type flow is evident by the 
good agreement between mean C P distributions for a single cylinder from Phases 1 and 2 and data from Roshko 19 
acquired at a Reynolds Number of 8.6 x 10 6 . 

C. Measurement Techniques 

The measurement techniques used in this study were selected based on their ability to meet the primary objective 
of each phase. In Phase 1, 2-D Particle Image Velocimetry provided detailed insight into the steady and unsteady 
near-field flow structure. The process used to acquire and process the PIV data is described in Ref. 12. Since the 
objective of Phase 2 was to document and characterize the flow on the surface of each cylinder, steady and unsteady 
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pressures were measured on the front and rear cylinders simultaneously using the static orifices and transducers 
described in the previous section. An electronic scanning pressure system was used to acquire the steady pressures, 
which were then converted to pressure coefficients. The uncertainty in C P is estimated to be ±0.02. For the unsteady 
pressures, signals from the piezoelectric resistive transducers were amplified and then sent to the A/D converter 
where they were AC-coupled, low-pass filtered at 10 kHz, and digitized using a sample rate of 25.6 kHz and sample 
duration of 32 seconds. To check the calibration coefficients and monitor the health of the transducers during the 
test, the signals for the transducers in the circumferential arrays were split before the A/D converter so the total 
signal could be digitized at a slower rate, averaged, and compared to the static pressure data. In general, this check 
showed that the transducers were quite stable and experienced little drift during the course of the test. 

To help interpret the unsteady pressure data, flow visualization was performed to identify surface flow features. 
These features were visualized by applying a mixture of titanium dioxide (Ti0 2 ) and kerosene to the model surface. 
As the tunnel was started and adjusted to the proper operating conditions, the mixture was transported along surface 
streaklines. Over time, the kerosene evaporated and the Ti0 2 was left on the surface in a solidified form, revealing 
the surface flow pattern. 


III. Discussion of Results 


A. Consistency Between Phases 

One major concern in conducting multi-phase wind tunnel experiments is insuring that tunnel conditions and the 
physical features of the flow remain consistent between tunnel entries. Consistency is very important in this case 
since we desire to combine off-body information obtained via PIV from Phase 1 and surface flow information from 
Phase 2 to develop a comprehensive picture of the interference between tandem cylinders. There are some notable 
differences between the two phases of this study. Firstly, the cylinder diameter was increased in Phase 2 to 
accommodate the unsteady pressure transducers and secondly, the tunnel velocity was decreased to maintain the 
same Reynold’s number. To assess consistency, C P distributions acquired on a single cylinder during each phase 
were compared. Figure 3, which was used to show the effectiveness of the transition strip, also shows the good 
agreement between the pressures acquired during the two phases. The only notable differences occur near the edges 
of the transition strips whose location is indicated by the two, gray vertical columns. To identify global differences 
in the flow between the two phases, the drag coefficient was computed by integrating the pressures measured on 
each cylinder between the limits 0 = 0 deg. and 0 = 360 deg. Figure 4 shows a comparison of the drag coefficients 
for each L/D tested. The overall trends are consistent; however, the drag on the front cylinder is slightly higher in 
Phase 2. This is also the case for the rear cylinder but the differences occur mainly at L/D’s between 3 and 4. 
Despite these differences, the consistency was deemed adequate to characterize the flow at L/D = 1.435 and L/D = 
3.7 using data from both phases. 

B. Flow Characterization 

Although several L/D spacings were tested, this discussion will focus on characterizing the flow associated with 
L/D = 1.435 and L/D = 3.7. For comparison, results for a single cylinder will also be presented. The flow field 
description will refer to the gap region and the near-wake region. The gap region includes flow on the rearward face 
of the front cylinder, on the forward face of the rear cylinder, and between the two cylinders. The near-wake region 
includes flow on the rearward face of the rear cylinder and the immediate vicinity downstream of the rear cylinder. 
For the PIV data presented, distances have been normalized by the cylinder diameter and are referenced to the center 
of the front cylinder for the gap flow figures and to the center of the rear cylinder for figures featuring the near-wake 
flow. Color contours of the spanwise vorticity have been normalized by the freestream velocity and cylinder 
diameter. 

1. Gap Region - Mean Flow 

The flow in the gap region is of primary importance because its interaction with the rear cylinder greatly 
influences the drag and the structure of the near- wake flow. Figures 5 and 6 show the mean flow structure in the gap 
region for L/D = 1.435 and 3.7 using streamlines computed from the streamwise and vertical velocity components. 
In the figures, the streamlines have been superimposed over color contours of the streamwise velocity. For L/D = 
1.435, the flow separates from the front cylinder and impinges on the rear cylinder to form a single cell of 
recirculating flow in the lower portion of the gap. Computations by Khorrami, et al. 13 have shown that this type of 
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flow pattern can result from flow angularity as small as 0.5 degrees. An extensive effort was made during both 
phases to align the model to the flow using both single and tandem cylinder configurations. The process was 
interative and involved adjusting the sideslip angle of the model in 0. 1 deg. increments while checking the symmetry 
between pressures on the upper and lower surfaces of the cylinder. Even after finding the position that produced 
good symmetry for the single cylinder and L/D = 3.7, the asymmetric flow pattern still persisted at L/D = 1.435. 
The existence of this flow pattern in both phases of this study as well as in other experiments at both subcritical 5 6 and 
supercritical 7 Reynolds numbers suggests that either there may be other contributing factors to the asymmetry 
besides sideslip angle or this may be a preferred flow state. In contrast, the mean flow pattern for L/D = 3.7 is 
typical for a cylinder in a crossflow; namely, a region of recirculating flow behind the cylinder comprised of two, 
counter-rotating cells. For these particular cylinder spacings, the differences in the gap flow have a greater effect on 
the C P distribution for the rear cylinder than for the front cylinder. This is evident in Figs. 7 and 8 which show mean 
C P distributions for the front and rear cylinders. Despite slight differences in the pressure recovery region of the 
front cylinder, the C P distributions for L/D = 1.435 and L/D = 3.7 are almost identical. They both show a higher 
base pressure than the single cylinder and both show separation at 0= 110 deg. and 6= 250 deg. The C P 
distributions shown in Fig. 8 for the rear cylinder have radically different trends on the forward face between 0 = 0 
deg. and 0 = 90 deg. and 0 = 270 deg. and 6= 360 deg. For L/D = 1.435, the flow is dominated by the attachment 
of the shear layer and the associated flow bifurcation. For L/D = 3.7, the flow stagnates at 0 = 0 deg. so the trend is 
similar to the front cylinder although the pressure level is lower. Differences are also noted on the rearward face of 
the cylinder between 0 = 130 deg. and 0 = 230 deg. Here, the pressure recovery is greater for L/D = 1 .435 than L/D 
= 3.7, which produces the negative drag coefficient shown in Fig. 4. 

In addition to streamlines and surface pressures, surface flow patterns shown in Figs. 9 and 10 provide further 
insight into the interaction of the gap flow with the cylinders. The photographs were taken with the camera 
perpendicular to the X-Z plane and pointing in the negative Y direction. As such, the orientation of the cylinders 
differs from that shown in Figs. 5 and 6. The front cylinder is on the right, the rear cylinder is on the left, and the 
flow is from right to left. For L/D = 1.435, the surface flow pattern on the front cylinder in Fig. 9 shows the 
transition strip (1) and streaks from the streamwise vortices (2) it generates. These vortices coalesce along the 
primary separation line (3). It also shows flow moving in the spanwise direction (4) between the primary separation 
line (3) and secondary separation line (5). Since the shear layers bound the flow, spanwise flow is necessary to 
maintain the momentum balance and may explain the streamline pattern near the upper shear layer in Fig. 5 where 
the recirculating flow traveling toward the front cylinder appears to change direction abruptly near X/D = 0.7 and 
Y/D = 0.5. For L/D = 3.7, a similar surface flow pattern is seen on the front cylinder in Fig. 10 however there is 
little evidence of spanwise flow between the primary separation line (3) and secondary separation (5) line. If 
spanwise flow is characteristic of tandem configurations with small spacings, then computer simulations of such 
flows will have to be 3-D and have some estimate of the spanwise correlation length to predict these flows 
accurately. 

For the rear cylinder, the surface flow patterns in Figs. 9 and 10 show little difference. The primary separation line 
(6) appears to occur at the same location for both configurations but shows more uniformity along the span of the 
cylinder for L/D = 3.7. From the C P distributions shown in Fig. 8, the flow appears to separate near 0 = 120 deg. 
and 6 = 240 deg for L/D = 1.435 and near 6 = 130 deg. and 6 = 230 deg for L/D = 3.7. This is aft of the separation 
location estimated for the front cylinder. The difference in the separation location between the front and rear 
cylinders may be due to the differences in the turbulence levels in the boundary layer. The front cylinder, boundary 
layer undergoes transition but the rear cylinder, immersed in an unsteady flowfield, has a boundary layer that is 
already turbulent due to turbulent flow attaching to the forward surface. 

2. Gap Region - Unsteady Flow 

Full characterization of the unsteady flow in the gap region requires information about the instantaneous flow 
structure and how this structure interacts with the surface of the cylinders. Figures 1 1 and 12 compare the fluctuating 
pressure coefficients on a single cylinder to the two tandem configurations. The fluctuating pressure coefficient was 
computed by integrating the spectra across the frequency range to obtain the mean square value, taking the square 
root, and then normalizing by the freestream dynamic pressure. The relative uncertainty in the computed spectra 
and fluctuating pressure coefficient are estimated to be 10% and 5.4%, respectively. Similar to the mean Cp 
distributions, Fig. 11 shows some variation between the three cases but in general, the unsteady gap flow appears to 
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have minimal impact on the pressure fluctuations on the front cylinder at these conditions. The only effects are seen 
in the vicinity of 6 = 90° and 6 = 270° due to unsteadiness in the flow separation point and in the vicinity of 0 = 
135° and 6 = 225° where vortices formed behind the cylinder interact with the surface. Sun, et al. 8 performed 
similar measurements on tandem cylinders at a supercritical Reynolds number of 6.5 x 10 5 and found that the 
fluctuating pressures on the front cylinder collapsed into a single curve for L/D < 4. Although the curves in this 
study did not collapse, they do show similar trends and only slight variations in pressure level. In addition to the 
pressure distributions in Fig. 3, this may provide further evidence of the transition strip’s effectiveness in producing 
post-critical flow on the cylinders. 

Compared to the front cylinder, the effects of unsteady flow in the gap region on the rear cylinder are more 
pronounced. Overall, the fluctuating pressure levels in Fig. 12 are higher for L/D = 3.7 than for L/D = 1.435. For 
the former case, the maximum value of C P r measured using the circumferential tranducers is 0.68 and occurs in the 
vicinity of 6= 45° and 6= 315°. The maximum levels for L/D = 1.435 are 0.14 and 0.18 and occur in the vicinity 
of 0 = 60° and 0 = 300° where the shear layers attach to the rear cylinder. These surface flow characteristics are 
corroborated by representative frames of the instantaneous spanwise vorticity shown in Figs. 13 and 14. Figure 13 
shows the attachment of the lower shear layer on the rear cylinder for L/D = 1.435 and Fig. 14 shows the formation 
of a vortex and the interaction of the vortex sheet with the rear cylinder for L/D = 3.7. In both cases, smaller 
vortical structures formed in the shear layer can be seen in proximity of the rear cylinder, which could possibly 
contribute to higher surface pressure fluctuations and broader peaks in the spectra. 

Spectra for individual transducers on the front and rear cylinders at L/D = 1.435 and L/D = 3.7 were examined to 
identify frequencies associated with features such as flow separation, shear layer attachment, and vortex shedding. 
Figures 15 and 16 show pressure power spectra computed for transducers around the circumference of the front and 
rear cylinders. Since the spectra are fairly symmetric between the upper and lower halves of each cylinder, spectra at 
0= 90, 135, and 180 deg. on the front cylinder and 0= 0, 45, and 90 deg. on the rear cylinder have been compared 
in Figs. 17-20 for this discussion. 

On the front cylinder, the spectra for L/D = 1.435 in Fig. 17 and the spectra for L/D = 3.7 in Fig. 18 have 
comparable energy levels and features with the exception of the peaks in Fig. 18 associated with vortex shedding. At 
6 = 90 deg., the pressure level increases between 1 kHz and 10 kHz. After examining the surface flow patterns, this 
increase is believed to be due to the streamwise vortices generated by the transition strip as identified in Fig. 9. In 
comparison, the spectra on the rear cylinder in Figs. 19 and 20 at 6 = 90 deg. do not show this feature. At 6 = 135 
deg., the spectra for both configurations show a broadband “hump” centered at 1 kHz. Although not shown in this 
paper, the spectra for the transducers in the spanwise array at this same azimuth all have this feature as well. It was 
also seen in spectra from a hot-wire survey conducted through the shear layer at X/D = 0.5 during Phase 1 and 
appears in the spectrum on the rear cylinder at 6 = 45 deg., which is near the location where the flow attaches. This 
provides conclusive evidence that this hump is associated with the shear layer. Based on an examination of 
instantaneous PIV images, it is unclear whether the hump is reflecting the movement of the shear layer, the vortical 
structures of various scales being convected, or both. In addition to the “hump”, Fig. 18 also shows peaks at 190.6 
Hz and 381 Hz associated with vortex shedding from the front cylinder. For both configurations, the spectra at 6 = 
1 80 deg. are rather benign. 

Spectra in Figs. 19 and 20 for the rear cylinder show similar features but the pressure levels are much higher 
than those measured on the front cylinder over the entire frequency range. For L/D = 1 .435, the “humps” are present 
at both 6= 45° and 6= 90° but the energy level is higher at 6= 45° due to its proximity to the attachment point. 
For L/D = 3.7, peaks associated with vortex shedding as well as harmonics of this frequency show up in the each of 
the sensors. At 6= 0°, more energy is in the harmonic peak than in the primary shedding peak. Since this 
transducer is located at the centerline between the front and rear cylinders, it senses the vortices shed from both 
sides of the front cylinder instead of just one. Similar features were seen in hot-wire spectra from traverses behind 
the rear cylinder. 
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3. Near-Wake Region 

Although the near-wake behind the rear cylinder is assumed to be dominated by vortex shedding, the spacing 
between the cylinders and consequently the gap flow affects the mean and unsteady wake structure as well. This is 
evident in Figs. 21a-c and Figs. 22a-c which show mean streamlines and a representative frame of the instantaneous 
spanwise vorticity behind a single cylinder and the rear cylinder of the two tandem configurations. In Fig. 21 the 
size of the recirculation bubble in terms of width and length differs between the three configurations due to 
differences in separation location discussed earlier. Additional flow field measurements at other L/D’s are needed to 
determine if the wake width and closure length are constant for a given flow regime. The size of this region has 
implications on pressure recovery and indicates where the unsteady flow is most likely to interact with the cylinder 
surface. 

In terms of the unsteady wake structure, Figs. 22a-c show that the wake is more organized for the single cylinder 
and L/D = 3.7 than for L/D = 1.435. From the mean streamlines and steady surface pressures, it is apparent that for 
L/D = 1.435 flow acceleration around the rear cylinder is suppressed, and consequently, pressure recovery is 
reduced and base pressure is higher. This, in turn, may lead to a less organized vortex street in the wake, as both the 
positive and negative levels of vorticity are reduced. Lower vorticity levels would explain why vortex shedding 
peaks are absent from the spectra in Fig. 19 for L/D = 1.435. To investigate this further, pressure spectra acquired 
on the rear cylinder at 6 = 225 deg. and velocity spectra acquired during Phase 1 downstream of the rear cylinder 
were plotted versus Strouhal number. Figure 23 shows the measurement locations relative to the rear cylinder and 
Figs. 24 and 25 show pressure and velocity spectra. Both measurement techniques capture the shedding from the 
single cylinder and L/D = 3.7 but for L/D = 1.435, no narrowband peaks appear in the pressure spectra. A broadband 
peak is distinguishable in the velocity spectra at St = 0.280 but it is unclear from the experimental data whether this 
is due to a less-organized vortex sheet or due to “undulations” in the shear layer. An examination of additional 
frames of instantaneous spanwise vorticity did not show evidence of vortex formation behind the rear cylinder for 
L/D = 1.435 but did show the type of undulations seen in free shear layers. These results indicate that the 
interference effects for the L/D = 1.435 configuration, and perhaps other configurations with small spacings, can be 
considerable and the assumption that it sheds like a single cylinder may be invalid. 

In addition to spectra, the pressure correlation length is an important parameter for characterizing unsteady flow 
on cylinders. Not only does it indicate the distance over which the surface pressure fluctuations are correlated, but it 
also provides information about vortex shedding. In research on single cylinders, Bruun, et al. 20 found that longer 
correlation lengths (i.e., distances over which correlation coefficients are significant) are indicative of coherent 
vortex shedding continuing over longer distances downstream from the generating cylinder. Conversely, shorter 
correlation lengths indicate smaller distances over which the wake vortices stay coherent. 

The broadband correlation length was computed on the single cylinder and rear cylinder of the tandem 
configurations using the spanwise row of transducers at 0-135°. Starting with the transducer closest to the floor of 
the tunnel as the reference point, the correlation coefficient was computed between it and each transducer above it. 
This process was repeated using each successive transducer as the reference point until correlations had been 
computed for every possible transducer combination. In Figs. 26 and 27, the correlation coefficients computed using 
different reference points along the span of the cylinder have been plotted for L/D = 1.435 and L/D = 3.7. The 
correlation length for L/D = 1.435 is much shorter than for L/D = 3.7. This is indicative of a wake that shows 
coherent vortex shedding over a very short distance. The correlation length for L/D=3.7, on the other hand, is much 
longer, indicating coherent vortex shedding over longer distances downstream. This is confirmed by the PIV images 
shown in Fig. 22. The two plots also provide insight into the dependence of the correlation, and consequently 
correlation length, on the location of the reference transducer. For L/D = 1.435, the curves collapse indicating that 
the correlation is independent of the location of the reference point. In contrast, the curves for L/D = 3.7 begin to 
diverge after Z/D = 1.25 which shows the dependence of correlation on the location of the reference point and 
suggests that the flow may be nonuniform along the cylinder span. Howell and Novak 21 documented similar results 
for a single cylinder but also found that the correlation length was dependent on end effects. Although further 
analysis is needed, this may explain the differences in Fig. 27 given that the points where the curves differ most 
correspond with transducers located near the tunnel ceiling. 
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To examine the relationship between correlation lengths based on pressure fluctuations and correlation lengths 
based on velocity fluctuations, results from the current phase are compared to correlations from Phase 1 in Fig. 28. 
The fluctuating velocities were measured using two hot wires at the same location behind the rear cylinder where the 
velocity spectra were acquired. The first hot wire was fixed at mid-span while the second hot-wire was traversed 
toward the tunnel ceiling. For the fluctuating pressure correlations, the transducer at mid-span was used as the 
reference point. With the exception of L/D = 1.435, the velocity correlations in Fig. 28 are lower than the pressure 
correlations, at a given separation, indicating a shorter correlation length. This is to be expected due to vortex 
shedding and increased three-dimensionality of the flow as the wake develops. The differences between the 
correlations for L/D = 1.435 may be due to spanwise flow on the surface of the rear cylinder between the primary 
and secondary separation lines. 


Conclusion 

A multi-phase experimental study has been conducted to characterize the unsteady interaction between two 
cylinders in tandem arrangement. Transition strips applied to the upper and lower surface of the front cylinder were 
successful in achieving post-critical flow around the cylinder. Steady and unsteady surface pressures, flow 
visualization, PIV, and hot-wire anemometry have been used to establish a comprehensive dataset for the flow 
around a single cylinder and a tandem configuration for L/D = 1.435 and L/D = 3.70. Characteristics about the flow 
in the gap and near- wake regions are summarized below. 

1. Mean flow in the gap region at L/D = 1.435 exhibited a single recirculation cell and impingement of the free 
shear layer on the front of the rear cylinder. Significant spanwise flow appeared to occur between the primary and 
secondary separation locations. At L/D = 3.7, two counter-rotating recirculation cells were generated, and 
stagnation occurred on the front face of the rear cylinder at 0 = 0 deg. There was little evidence of spanwise flow 
between the primary and secondary separation locations. Separation on the rear cylinder appeared to occur about 10 
degrees further forward for L/D=1.435 than for L/D=3.7. 

2. Unsteady flow in the gap region was characterized by much higher unsteady pressures on the forward face of 
the rear cylinder than on the rearward face of the front cylinder, due, most likely, to the impingement of a turbulent 
shear layer on the rear cylinder. Power spectra from transducers on the surface of the front cylinder showed peaks at 
vortex shedding frequencies for the single cylinder and L/D=3.7. Vortex shedding for L/D=1.435 was suppressed. 
There was also evidence of increased spectra levels at high frequencies, at 6 = 90 and 6 =210 deg, possibly due to 
the streamwise vorticies shed by the transition strips, as seen in the surface flow visualization. Spectra levels on the 
forward face of the rear cylinder were 3 to 5 orders of magnitude greater than on the front cylinder, likely due to the 
impingement of the turbulent shear layer separated from the front cylinder. 

3. Flow in the near- wake region was characterized by vortex shedding. However, the wakes for the single 
cylinder and rear cylinder at L/D=3.7 were more organized than for L/D=1.435. Acceleration around the rear 
cylinder for this L/D was suppressed, pressure recovery was reduced, and the base pressure was higher, which may 
have contributed to the lack of organization in the wake. No evidence of a vortex shedding peak in the spectra for 
this L/D existed. The existence of a slight “hump” in the spectra measured by a hot wire in the wake is surmised to 
be due to an “undulation” of the shear layers rather than a real vortical structure. Correlation lengths on the rear 
cylinder for L/D = 3.7 are approximately equal to the single cylinder whereas the correlation length for L/D = 1.435 
is shorter. Of the two tandem configurations discussed, the near-wake structure for L/D = 3.7 is most like that of a 
single cylinder. 

These results, along with those obtained during the other phases of this study, will be used to establish a 
comprehensive dataset consisting of surface and off-body measurements on a tandem cylinder configuration. Such 
data can be used to 1) characterize the interaction of the wake of the front cylinder with the rear cylinder for 
different spacings, 2) examine differences between the wake of the entire tandem configuration and a single cylinder 
due to this interaction and 3) develop and validate computational techniques for the prediction of unsteady flows. 
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Figures 



Figure 1. Tandem cylinder configuration in the 
BART facility. 



Figure 2. Model schematic and coordinate 
reference system. 
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Figure 3. Pressure distribution on a single cylinder with transition strip. 



Figure 4. Effect of cylinder spacing on drag coefficient. 
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Figure 6. Mean streamlines in gap region (L/D = 3.7). 



Figure 7. Mean pressures on front cylinder. 



Figure 8. Mean pressures on rear cylinder. 


11 

American Institute of Aeronautics and Astronautics 



Figure 9. Surface flow visualization (L/D = 1.435). 



Figure 10. Surface flow visualization (L/D = 3.7). 
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Figure 11. Fluctuating pressure coefficient on front Figure 12. Fluctuating pressure coefficient on rear 
cylinder. cylinder. 
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Figure 13. Instantaneous spanwise vorticity (L/D = 1.435). 
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Figure 14. Instantaneous spanwise vorticity (L/D = 3.70). 
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Figure 15. Pressure power spectra on front cylinder for L/D = 1.435 and L/D = 3.7 compared to single 
cylinder. 
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Figure 16. Pressure power spectra on rear cylinder for L/D = 1.435 and L/D = 3.7 compared to single 
cylinder. 
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Figure 17. Spectra on front cylinder (L/D = 1.435). Figure 18. Spectra on front cylinder (L/D = 3.7). 




Figure 20. Spectra on rear cylinder (L/D = 3.7). 
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Figure 21. Mean streamlines in near-wake region 


c) L/D = 3.70 

Figure 22. Instantaneous spanwise vorticity. 
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Figure 23. Pressure and velocity spectra measurement locations. 
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Figure 24. Pressure spectra on rear cylinder at 
0 = 225 deg. 



Figure 26. Fluctuating pressure correlation on rear 
cylinder (L/D = 1.435) 



Figure 25. Velocity spectra behind rear cylinder at 
X/D = 2.74, Y/D = -0.685. 



Figure 27. Fluctuating pressure correlation on rear 
cylinder (L/D = 3.7) 
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Figure 28. Correlations computed from hot-wire and surface pressure measurements. 
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